Spin-polarized Josephson current in S|F|S junctions with inhomogeneous magnetization 
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We study numerically the properties of spin- and charge-transport in a current-biased nanoscale diffusive 
superconductor] ferromagnet] superconductor junction when the magnetization texture is non-uniform. Specif- 
ically, we incorporate the presence of a Bloch/Neel domain walls and conical ferromagnetism, including the 
role of spin-active interfaces. The superconducting leads are assumed to be of the conventional s-wave type. 
In particular, we investigate how the 0-n transition is influenced by the inhomogeneous magnetization texture 
and focus on the particular case where the charge-current vanishes while the spin-current is non-zero. In the 
case of a Bloch/Neel domain-wall, the spin-current can be seen only for one component of the spin polarization, 
whereas in the case of conical ferromagnetism the spin-current has the three components. This is in contrast to 
a scenario with a homogeneous exchange field, where the spin-current vanishes completely. We explain all of 
these results in terms of the interplay between the triplet anomalous Green's function induced in the ferromag- 
netic region and the local direction of the magnetization vector in the ferromagnet. Interestingly, we find that 
the spin-current exhibits discontinuous jumps at the O-tt transition points of the critical charge-current. This 
is seen both in the presence of a domain wall and for conical ferromagnetism. We explain this result in terms 
of the different symmetry obeyed by the current-phase relation when comparing the charge- and spin-current. 
Specifically, we find that whereas the charge-current obeys the well-known relation Ici<t>) ~ — /c(27r — (p), the 
spin-current satisfies Is (<?!>) = /s(27r — (p), where is the superconducting phase difference. 

PACS numbers: 74.78.Na 



I. INTRODUCTION 

Because of the interesting phenomena that 
superconductor! ferromagnet I superconductor (S|F|S) struc- 
tures exhibit, including their poten tial a ppHcations in 
spintronics^ and quantum computing^'*'^'*', this field of 
research is presently studied extensively.*^ Usual electronic 
devices are based on the properties of flowing electrons 
through circuits, whereas spintronic devices are based on 
direction and number of flowing spins. In many spintronics 
devices, like magnetic tunnelling junctions, spin polarized 
currents are generated when an imbalance between spin 
up and down carriers occurs. This imbalance can arise 
e.g. by using magnetic materials or applying a magnetic 
field. The discovery of the giant magnetoresistance (GMR) 
effecP today forms the basis of the leading technology for 
information storage by magnetic disc drives. Spin coupling 
and its advantageous high speeds at very low powers^ of 
these devic es prom ise applications for logic and storage 
applications.Ii^lilllSI 

The possibility of a 7r-state in a S|F|S systems was pre- 
dicted theoretically in Refs.'^'"* and has been observed 
experimentally.'^ Near such a transition point, the junction 
ground state energy has two minima versus cf) at (j> = and 
(f> = TT. The coexistence of stable and metastable and tt 
states in the transition zone can produce two flux peaks for 
one external quantum flux in superconducting quantum inter- 
ference device (SQUID)-like geometry, and renders the sys- 
tem a qubit.'J^ The characteristic length of the ferromagnetic 
layer where the first O-vr transition occurs is of the order of 
the magnetic coherence length ^p. In the dirty limit, that is 
achievable in most of the experimentally studied S|F|S struc- 



tures. Here, S^p is given by jh where D denotes the dif- 
fusion constant and h is the magnitude of ferromagnetic ex- 
change field. Therefore, the experimental observation of such 
O-vr transitions in nanoscale devices requires a low exchange 
energy h. Such conditions were achieved using weak ferro- 
magnetic CuNi or PdNi alloys. The experimental observa- 
tions of the critical charge-current oscillations shows such 0- 
TT transition s as a fun ction of the ferromagnet thickness and 
temperature .EEEIEII The consequence of the exchange split- 
ting at the Fermi level^^is that the Cooper pairs wave function 
shows damped oscillations in the ferromagnet, resulting in the 
appearance of the well known 7r-state in S|F|S systems.''' In 
contrast to the usual 0-state in superconductor-normal metal- 
superconductor junctions, the phase shift equal to tt across 
the junction in the ground state reverses the direction of the 
supercurrent,tS]and considerably changes the density of states 
(DOS) in the F metal. '^ The 7r-states can also be observed in 
nonmagnetic junctions of high-Tc superconductor^^ and in 
non-equilibrium nanoscale superconducting structures.*^ 

In the ballistic limit, the transport properties of a S|F|S 
junction can be understood on a microscopic level in terms 
of Andreev bound-states.^^ The O-tt transition is then due to 
the spin dependence of the Andreev bound states.l^S Because 
of the averaging of the quasiclassical Green's function^^ over 
momentum directions, the relevant equations simplify in the 
dirty transport regime. This averaging of Green's function 
can be understood by noting that in the presence of impuri- 
ties and scattering centers, the direction of motion of elec- 
trons are random and physical quantities should be averaged 
over all directions. This averaging is valid as long as the 
mean free path of the diffusive layer is much smaller than 
length scales of the system that are superconducting coher- 
ence length = -\/_D/Ao and the decay length of Cooper 
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pairs wave functions inside ferromagnet = \/D/h, where 
Aq is the superconducting order parameter. The charge- 
current Ic{4') ™d the local DOS are two principle quantities 
that are strongly influenced by the proximity effect. These 
two quantities were studied for various geometries by using 
of quasiclassical Green's function in clean and dirty limits in 

several works, C. g P6I27I28I29I30I3 Il32l33l34l35l36l^ 

Up to now, the majority of works studying S|F|S junctions 
have considered a homogeneous exchange fie ld in the fe r- 
romagnet, including half-metallic ferromagnets.'^^ElIlIEoEIlin 
the presence of inhomogeneous magnetization textures, sev- 
eral new effects have been predicted in the literature includ- 
ing the possibility of a long-range triplet component. Such 
an inhomogeneous magnetization texture may be created ar- 
tificially by setting up several layers of ferromagnets with 
misaligned magnetizations."*^ Alternatively, inho- 
mogeneous magnetization may arise naturally in the presence 
of domain walls or non-trivial patterns for the local ferro- 
magnetic moment. An example of the latter is the conical 
ferromagnet Ho. Very recently, two theoretical studies have 
predicted qualitatively new effects in S|F and S|F|S hybrid 
structures where F is a conical ferromagnet"*'^^''. Due to the 
inhomogeneous nature of the magnetization in Ho, the spin- 
properties of the proximity-induced superconducting correla- 
tions are expected to undergo a qualitative change compared 
to the case of homogeneous ferromagnetism. Such changes 
may also be expected in the domain-wall case. A more re- 
alistic modeling of hybrid structures involving superconduc- 
tors and ferromagnets demands that such non-trivial magneti- 
zation textures and a lso the spin-dependent properties of the 
interface region^^^El are taken into account seriously. It was 
recently shown that the latter may induce qualitatively new 
features in the local DOS of S|F layers^^ and S|N layers with 
magnetic interfaceJ^. 

Another consequence of inhomogeneous magnetization, be 
it in the form of multiple misaligned layers or intrinsic non- 
uniformity within a single ferromagnetic layer, is that the 
Josephson current should become spin-polarized. This has 
been noted by several authors in the context of superconduc- 
tors coexisting with helimagnetic or spiral magnetic ordeiE^^^ 
as well as ferromagnetic superconductorsSIS] However, the 
spin-polarization of the Josephson current has not been stud- 
ied in the arguably simplest case of a single ferromagnetic 
layer with inhomogeneous magnetization contacted by two 
conventional s-wave superconductors. 

To this end, we develop in this paper a model for an S|F|S 
junction where both inhomogeneous magnetization and spin- 
active interfaces are incorporated, and then proceed to solve 
the problem numerically. More specifically, we will investi- 
gate variations of spin- and charge- currents versus changing 
of the thickness of F layer dp for a hybrid S |F| S structure with 
s-wave superconductors. We find that a spin-current flows 
through the junction whenever the magnetization is inhomo- 
geneous, and that it features discontinuous jumps whenever 
the junction undergoes a O-n transition. We compare these 
variations for three types of magnetization textures i.e., homo- 
geneous, domain wall, and a conical exchange field. We also 
show that for certain values of dp, the critical charge-current 



vanishes whereas a pure spin-current flows through the sys- 
tem. Moreover, we demonstrate how it is possible to obtain a 
pure spin-current by tuning the phase difference between the 
s uperconductors . 



II. THEORY 

To investigate the behavior of the ferromagnetic Josephson 
junction, we employ a full numerical solution of the quasiclas- 
sical equations of superconductivity^^ in the diffusive limit^^, 
which allows us to access the full proximity effect^^ regime. 
Importantly, we also take into account the spin-dependent 
phase-shifts (spin-DIPS) microscopically that are present at 
the superconductor I ferromagnet interfaces. For the purpose of 
stable and efficient numerical calculations, it is convenient to 
employ the Ricatti-parametrization of the Green's function as 
followsPJ^ 



'Af{i - 77) 27y;7 



Here, = I since 



AC= (1 + 77)"^ 7v:= (1 + 77)" 



(1) 



(2) 



We use for 2x2 matrices and . ? . for 4x4 matrices. In 
order to calculate the Green's function g, we need to solve 
the Usadel equatiorP^ with appropriate boundary conditions 
at X ~ —dp/2 and x — dp/2. We introduce the supercon- 
ducting coherence length as ^5 = y^Ds/Ao- Following the 
notation of Ref.l^, the Usadel equation reads 

Dd{gdg) + i[ep3 + dmg[h ■ai,{h- af],g] - 0, (3) 

and we employ the following realistic boundary conditions for 
all our computations in this paperP^l 

2CdFgdg = [sbcs (</>),. 9] + iIGs/Gt) [diag(T3, 73), g] (4) 

at X = —dp/2. Here, d = and we defined C — Rb/Rf as 
the ratio between the resistance of the barrier region and the 
resistance in the ferromagnetic film. The barrier conductance 
is given by Gt, whereas the parameter G5 describes the spin- 
DIPS taking place at the F side of the interface where the mag- 
netization is assumed to parallel to the z-axis. The boundary 
condition at x — dp/ 2 is obtained by letting G5 (^^5) 
and gBCs(0) ^ [-^bcs in Eq. where 



7Bcs('^)=lI2s/(l + c)e'^/^ 



(5) 



Above, Gg is allowed to be different from G5 in general. 
For instance, if the exchange field has opposite direction at 
the two interfaces due to the presence of a domain wall, one 
finds Gs = -Gs- The total superconducting phase differ- 
ence is (f), and we have defined s — sinh(z9),c = cosh(i?) 
with -d = atanh(Ao/£) using Aq as the superconducting gap. 
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Note that we use the bulk solution in the superconducting re- 
gion, which is a good approximation when assuming that the 
superconducting region is much less disordered than the ferro- 
magnet and when the interface transparency is small, as con- 
sidered here. We use units such that h — fc^ — 1. 

The values of Gs and Gt may be calculated explicitly from 
a microscopic model, which allows one to characterize the 
transmission {i^ ^} and reflection amplitudes {r^ cr} on the 
j E {S, F} side. Under the assumption of tunnel contacts 

and a weak ferromagnet, one obtains with a Dirac-like barrier 
modeP'5160 



Gt — Gr 



Gq — 2Gr 



and 



(6) 



(7) 



upon defining r„ = J2a Kir 

= Im{<^(r^J*}, r„^ = Im{<^«,)*}. 

For simplicity, we assume that the interface is characterized 
by N identical scattering channels. Omitting the subscript 'n', 
the scattering coefficients are obtained as: 



(8) 



with the definitions — + fc 

\/2m^Jj^', k^ = y/2mF{p,F + crh). (9) 



' ik^ Za and 



k' 



Here, = Zq + aZg is the spin-dependent barrier potential. 
Defining the polarization P — hj [ip vs\ the ferromagnet and 
the polarization ly = Zg /Zq for the barrier, we will set P = v. 

In this paper we will consider three types of inhomogeneous 
magnetic textures: Bloch, Neel and a conical structure. These 
structures are all different from a homogenous magnetic tex- 
ture. The first two types of magnetic textures are assumed 
to be located at the center of the F layer. The Bloch model 
is demonstrated by h = h{cos9y + sinOz) and its struc- 
ture is shown in Fig. [T] Similarly, the Neel model reads 
h = h{cos 9x + sin 9z) where we defined 6 as follows^-': 



— arctan(a;/c?vi/). 



(10) 



Here, dw is the width of domain wall and we assumed that the 
center of F layer is located at the origin, i.e, x — as shown 
in Fig. [1] 

For the conical case, we adopt a model where the magnetic 
moment rotates on the surface of a cone with defined apex 
angle a and turning angle /3. This structure is shown in Fig.l 
(a and /3 will determine the kind of material in use). If we 
assume that the distances of interatomic layers are a,^"* the 
spiral variation in the exchange field can be written as 

h = h{cosax + sma[sm{f3x/a)y + cos{(3x / a) z]) . (11) 

To characterize the transport properties of the system, we 
define the normalized charge and spin current densities ac- 
cording to: 



Ic.O 



(12) 
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FIG. 1: (Color online) i) The S|F|S junction with Bloch domain 
wall of ferromagnet and /() with conical type of ferromagnet. The 
magnetization texture for the Neel wall is obtained by replacing the 
a;-component of the magnetization with an y-component in case i). 
The blue arrows show the magnetic moments in F layer. The mag- 
netic moment for Bloch/Neel domain wall has two components and 
for conical type has three components. 



and 



r 

ls,0 



dsTr{p3f3(4)"}, (13) 



respectively, where e = s/Aq, x = x/dp, and — 
diag(l, 1, — 1, — 1). Here 1^ and P- are the charge- and the 
z-component of the spin-current flowing in the x-direction, 
respectively. The normalization constants are: 
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where Nq is the normal state DOS per spin and D is the dif- 
fusion constant. In general, the spin-current for other compo- 
nents of spin polarization j E {x, y, z} is given as: 



P 



deTr 







(15) 



Above, Pi, fi, and Tj are Pauli matrices that are defined 
in the appendix C and the reader may consult Appendix B 
for the derivation of the expression for II/I^ q. Under the 
assumption of an equilibrium situation, the Keldysh block of 
Green's function reads: 



r = [.r-.ntanh(/3£/2), 



(16) 



where g^and = —{psg^ps) are Retarded and Advanced 
blocks of g respectively, and (3 — 1/T is inverse temperature. 



III. RESULTS AND DISCUSSION 

We now present our main results of this paper, namely a 
study of how the critical currents depend on the thickness dp 
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of the junction in the presence of homogeneous and inhomo- 
geneous exchange field and spin-active interfaces. In order to 
focus on a realistic experimental setup, we choose the junction 
parameters as follows. For a weak, diffusive ferromagnetic al- 
loy such as Pdj;Nii_:r, the exchange field /i/Aq is tunable by 
means of the doping level x to take values in the range meV 
to tens of meV. Here, we will fix /i/Aq = 15, which typi- 
cally places the exchange field h in the range 15-25 meV. The 
thickness dp the junction is allowed to vary in the range 
dp/is e [0.5, 1.2], which is equivalent to 9 — 21.6 nm for 
a superconducting coherence length of 1^5 = 18 nm as can 
be obtained for e.g. Nb. This range of layer thicknesses dp 
are experimentally feasible.l^ The ratio of Gs/Gt is calcu- 
lated according to the microscopic expressions given in the 
previous section only for uniform and domain wall exchange 
fields because we will set G5 = for the conical ferromag- 
net. We choose = 1 eV and = 10 eV for the Fermi 
level in the ferromagnet and superconductor, respectively, and 
consider a relatively low barrier transparency of Zq = 3. The 
electron mass mp and tos in both the F and S regions is taken 
to be the bare one (~ 0.5 MeV). Any change in effective mass 
translates into an effective barrier resistance due to the Fermi- 
wavevector mismatch, which thus is captured by the param- 
eter Za. The interface region is assumed to exhibit a much 
higher electrical resistance than in the bulk of the ferromagnet, 
and we set ( = Rg/Rp = 4. For more stability in our com- 
putations we used the Ricatti parametrization and also inserted 
a small imaginary part S = 5 x 10^"^ Aq in the quasiparticle 
energy e, effectively modeling inelastic scattering. A consid- 
erable amount of CPU-time was put into the calculations of 
the current, as we solved for a fine mesh of both quasiparticle 
energies e and phase differences (f> for each value of the width 
dp. As will be discussed in detail below, we find that for S|F|S 
structures with spin-singlet s-wave superconducting leads, a 
spin-current exists only for domain wall structures and con- 
ical type of the ferromagnet layer, whereas it vanishes com- 
pletely in the case of a homogeneous exchange field. Both the 
charge- and spin-current are evaluated in the middle of the F 
region, a; = 0. The charge-current is conserved throughout the 
system, and its magnitude is thus independent of x. The spin- 
current, on the other hand, is not conserved and in fact suffers 
a depletion close to the S |F interfaces and vanishes completely 
in the superconducting regions. The critical charge-current is 
given by Ice —m3Xcf,{Ic{(f>)}, and the phase giving the critical 
current may be denoted We define the critical spin-current 
as Ics=Is{<t>c), which means that we are effectively consider- 
ing the spin-polarization of the critical charge-current, which 
should be the most sensible choice physically in a current- 
biased scenario. Note that this is different from the maximum 
value of the spin-current as a function of 



A. Critical currents vs. thiclcness dp for homogeneous 
exchange field 

First, we consider how the charge- and spin-currents are 
influenced by changing the thickness of F layer dp m the ho- 
mogeneous magnetic texture case. We fix the temperature at 
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FIG. 2: (Color online) The variation of the normalized critical 
charge-current versus the thickness of a homogeneous F layer. The 
inset panel zooms in on the behavior near the O-tt transition. As long 
as the exchange field is constant, we find that the spin-current Is 
vanishes. 



T/Tc ~ 0.2, and use the microscopic expression for spin- 
DIPS G^ at the two boundaries. The result is shown in the Fig. 
|2] The critical charge-current in the region of dp from 0.5^5 
to 1.2^5 vanishes at one point. This point is the first O-tt tran- 
sition point. We found that, for all strengths of the exchange 
field and spin-DIPS, the spin-current Is is zero. Unlike the 
case of spin-triplet superconductors, we can not see a ny sp in- 
current even for x and y directions of spin polarizatiorPSI^. In 
fact, one can confirm this finding analytically for all compo- 
nents of spin polarization at least for linearized Usadel equa- 
tion and transparent boundaries. The reason for the vanish- 
ing spin-current will become clear from the discussion in the 
following section, when noting that only the 5*^ = odd- 
frequency triplet and even-frequency singlet components are 
induced by the proximity effect in the ferromagnetic region. 



B. Critical currents vs. thickness dp for Bloch and Neel 
domain walls 

We now turn our attention to the first example of a non- 
trivial magnetization texture in the ferromagnet, namely the 
scenario of a Bloch or Neel domain wall. We use the same 
values for h and T as in the previous section, and set the do- 
main wall width dw to dw/dp = 0.5 and assume that it is 
centered in the ferromagnet. Although the domain-wall struc- 
ture dictates that the magnetization is not fully directed along 
the 2;-axis at the interfaces, we have verified numerically that 
the influence of the spin-DIPS parameter is negligible for 
our choice of parameters, such that we still can use the bound- 
ary conditions in Sec. |ll] 

The results of the variation of the normalized critical spin- 
and charge-currents vs. dp/£,s are shown in the Fig. [3], con- 
sidering here a Bloch wall texture. Contrary to the homoge- 
neous case considered in the previous section, we now see 
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FIG. 3: (Color online) The variations of normalized critical spin- 
and cliarge-currents vs. increasing the thiclmess dp of F layer with a 
Bloch domain wall structure. 



that a finite spin-current flows through the system. For this 
type of magnetization texture, we note that the spin-current 
only exists for one component of the spin polarization: the x- 
component. Only one component of the spin-current would 
be present also in the Neel domain wall case, as we shall ex- 
plain below. The spin-current features a discontinuous jump 
at the same value of the thickness where the charge-current 
undergoes a O-tt transition, namely dp/^s — 0.6. For this 
value of thickness the spin-current has a rapid variation. We 
have checked numerically with a very high resolution of dj? (a 
step of 5 X 10^'* for dp/^s) that this result does not pertain to 
noise or any error. For this type of magnetization texture, we 
also note that a spin-current only exists for one component of 
the spin polarization: the f -component. We will explain the 
reason for both the presence of such jumps in the spin-current 
and the polarization properties later 

We now explain why only one component of the spin- 
polarization is present both in the Bloch and Neel domain wall 
case. In order to understand the reason for this, it is instruc- 
tive to consider the interplay between the triplet anomalous 
Green's function /, given by 



fi - /t + /t) /n + /iT 



(17) 



and the local direction of the exchange field h. In S|F prox- 
imity structures, / tends to align as much as possible with h. 
For a homogeneous exchange field h in the z-direction, one 
thus obtains that only the 5*^ = opposite-spin pairing triplet 
component ft = is present, as is well-known. Consider 
now the Bloch domain wall case. The /-vector then contains 
only y- and z-components. Now, the spin expectation value 
of the Cooper pair is provided by 



(s)(xi(/xr), 



(18) 



the Neel domain wall case leads to the result that only a spin- 
polarization in the y-direction is present. Since we are eval- 
uating the spin-current in the middle of the F region, the z- 
component of the local exchange field is absent there. In 
that case, (S) should equal to zero according to our argument 
above. The reason for why a finite spin-current is neverthe- 
less obtained must be attributed to a lag between the / and h 
vectors, such that they do not follow each other exactly. One 
would expect that for a slower variation of the local exchange 
field, the lag would decrease. 



and we immediately infer that only a spin-polarization in the 
x-direction will be present. A similar line of reasoning for 



Critical currents vs. thickness dp for conical type of 
magnetization texture 



Finally, we turn our attention to the conical model for mag- 
netization, relevant to Ho. For simplicity, we set the Gs — 
at the two boundaries at —dp/^ and d^2. The distance be- 
tween the atomic layers a is equal to 0.02^;-, a = 47r/9, and 
rotating angle (3 — tt/6 per interatomic layer. These values of 
a, a and /3 are chosen based on the actual lattice parameters 
of Ho. The result of the investigation of how the critical spin- 
and charge-currents vary as a function of dp/^s is shown in 
the Fig. [4] In this case, we see a qualitatively new behavior 
for the charge-current as compared to the previous two sub- 
sections where we treated a homogeneous exchange field and 
a domain-wall ferromagnet, respectively. In Fig. |4] one ob- 
serves a superimposed pattern of fast oscillations on top of 
the usual O-tt oscillations, which are slower This is in agree- 
ment with the very recent work by Halasz et alM, who also 
reported the generation of rapid oscillations on top of the con- 
ventional O-TT transitions of the current in the weak-proximity 
effect regime. These faster oscillations pertain to the inho- 
mogeneous magnetization texture considered here, although 
they are not seen in the domain-wall case. This fact indicates 
that they are sensitive to the precise form of the magnetiza- 
tion structure in the ferromagnet, and that they do not appear 
simply as a result of a general inhomogeneity. 

As can be seen in the Fig. |4] the critical-charge current has 
five local minima, out of which three are O-vr transition points. 
In Fig. |4] the first dotted horizontal line indictates a minima 
which is irrelevant to a O-vr transition, whereas the three fol- 
lowing dotted lines indicate minima which correspond to such 
transitions. The last local minima is located near d^^/^s = 1.2 
and is not indicated by a dotted line in Fig. |4] This is in con- 
trast to the homogeneous and domain wall case, where only 
one O-TT transition point is seen in the range of dp considered 
here. As for the spin-current, the behavior is similar to the 
Bloch wall structure, with a rapid variation at the transition 
point. As mentioned previously, we have investigated these 
discontinuous jumps of the spin-current with a very high reso- 
lution for dp to ensure that do not stem from numerical errors 
or noise. We now proceed to an explanation for this effect. 
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FIG. 4: (Color online) Plot of the normalized critical spin- and 
charge-current vs. the normalized thickness dp of the F layer for 
a conical type of magnetization texture. All three components of the 
spin polarization have a considerable magnitude in the entire range 
of dp considered here. 

D. Origin of the discontinuous jumps in the spin-current 

In order to understand the mechanism behind the discon- 
tinuous jumps of the spin-current near the O-vr transition of 
the junction, we revert briefly to the original definition of the 
critical spin-current. It is defined as Is{4>c) where (pc is the 
value of the superconducting phase difference which gives the 
maximum (absolute) value of the charge-current. In effect, 
the critical spin-current is the spin-polarization of the critical 
charge-current, which is distinct from the maximum value of 
the spin-current. We now consider in detail the current-phase 
relation for both charge- and spin-transport near the transi- 
tion point located at dp/S^s — 0.772 (see Fig. The result 
for the current-phase relation is shown in Fig. |5l where we 
consider several values of dp near the transition point. From 
bottom to top, they range from dp/^g = 0.7655 to 0.7725 
in steps of 1 x 10"'^. A key point is that we have verified 
numerically that the charge-current is antisymmetric with re- 
spect to = TT whereas the spin-current is symmetric around 
this value. More specifically, whereas 

IM = -Wtt ~ (j>) (19) 

we find numerically that the spin-current satisfies 

=7,(2^-0) (20) 

This is consistent with the finding of Ref.l^ where transport 
between spin-triplet superconductors has been investigated. 
As a result, it suffices to restrict our attention to the range 
(j) £ [0,7r]. Next, we note that the charge-current is nearly 
sinusoidal to begin with (bottom curves of Fig. |5]l. Upon in- 
creasing dp, and thus approaching the transition point, higher 



FIG. 5: (Color online) Plot of the current-phase relation for the 
charge- and spin-currents. We are considering a conical magnetiza- 
tion texture, and the curves range from dp/^s = 0.7655 to 0.7725 
in steps of 1 x 10~^ in the charge-current panel. For the spin-current 
panels, the variation of the current-phase relation upon changing dp 
is negligible and we give results only for dp /(^s = 0.7725. Note that 
the critical spin-current nevertheless varies strongly with dp near the 
transition points as shown in Fig. |4] since the critical phase displays 
a strong dependence on ciF in this regime. 



harmonics in the current-phase relation become more protru- 
dent for the charge-current. However, the spin-current re- 
mains virtually unafffected by an increase in o?i?, and we plot 
the result only for dp/S^s — 0.7725. Upon increasing dp, the 
critical phase (pc moves away from tt/2 to lower values due 
to the presence of higher harmonics in the current-phase rela- 
tion. At the transition point, the phase jumps in a discontinu- 
ous manner to (pc > 7r/2 (dotted arrow in Fig. [S]). Now, the 
charge-current has a similar magnitude (in absolute-value) for 
this new value of cpc- The spin-current, on the other hand, has 
a different symmetry with respect to 4> as seen in Fig. |5] and 
varies less rapidly with dp. Therefore, the spin-polarization 
of the current makes a discontinuous jump at the transition 
point. 



IV. SUMMARY 

In summary, we have considered the transport of charge 
and spin in a nanoscale S|F|S Josephson junction when the 
magnetization texture is inhomogeneous in the ferromagnetic 
layer More specifically, we have investigated how charge 
and spin Josephson currents are affected by the presence of 
Bloch/Neel domain walls and conical ferromagnetism, includ- 
ing also the spin-active properties of the interfaces. We find 
that a spin-current flows through the junction whenever the 
magnetization is inhomogeneous, and that it features discon- 
tinuous jumps whenever the junction undergoes a O-n tran- 
sition. In the case of a Bloch/Neel domain-wall, the spin- 
current can be seen only for one component of the spin polar- 
ization (the component perpendicular to both the local direc- 
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tion of the exchange field and that of its derivative), whereas in 
the case of conical ferromagnetism the spin-current has three 
components. For a homogeneous exchange field, the spin- 
current vanishes. We explain the polarization properties of 
the spin-current by considering interplay between the triplet 
anomalous Green's functions induced in the ferromagnetic re- 
gion and the local direction of the magnetization vector in 
the ferromagnet. Moreover, we show how the discontinuous 
jumps in the spin-current stem from the different symmetries 
for the current-phase relation when comparing the charge- 
and spin-current. While the charge-current obeys the well- 
known relation /c(0) = — /c(27r — </)), the spin-current satis- 
fies Is{(j)) — Is{2TT—(f)), where is the superconducting phase 
difference. Our results open up new perspectives for appli- 
cations in superspintronics by exploiting Josephson junctions 
with non-homogeneous ferromagnets. 
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APPENDIX A: USEFUL RELATIONS FOR THE GREEN'S 
FUNCTION 

By introducing the auxiliary quantities 

da^Af = -Afmf, d^Af ^ -AfML, (Al) 
where we have defined 

^ - (9.7)7 + 7(5x7) , ^ = (9.7)7 + 7(9.7) , (A2) 

we find that the matrix derivative of the Green's function d^g^ 
has the following components: 

(9.ff^)ii = -VK - (1 - jjjAfVM, 
(9.ff^)i2 = 2(9.7)^:- 27^W, 

(9.3^)22 =M+(1-77)^W (A3) 

The indices above refer to particle-hole space, and each of the 
above elements is thus a 2 x 2 matrix in spin-space. 



APPENDIX B: QUASICLASSICAL EQUATION FOR THE 
SPIN-CURRENT 



mechanical expression for the expectation value of the spin- 
current: 

Usir)) = ^(Im{vI/t(r)V,diag(^,^*)vE'(r)}), (Bl) 
with a fermion operator basis 'i' given as 



vI/t(r) = (v4(r),V; (r),VT(r),^i(r)) 



(B2) 



Above, cr is the Pauli matrix vector. It should be noted that the 
spin-current jg is a tensor since it has a flow-direction in real 
space in addition to a polarization in spin-space. For clarity, 
we consider in what follows the (T2 -component corresponding 
to the polarization in the y-direction, as an example. We then 
get from Eq. (Bl i [using that Im{iz} = Re{z} for a complex 
number z] 

^■^'^^""^^ ^ ^Re{-(^|(r)V.^i(r)) + {^jlir)V (r)) 
+ (7/;T(r)V,Vl(r)) - (Vi(r-)VW4(r))} 



1 



4m 1 — >r' 
(^I(r)V'T(r')) + (^i(r)^|(r')) 

(V^t(^)W(^'))]- 



(B3) 



Using the notation of Ref.l^, we define the following repre- 
sentation for the Keldysh Green's function: 



(G^(r,r'; 



-iE(^3)™,([*(r-)„vI/t(r')„]- 

j 



(B4) 



It then follows from anticommutation that e.g. 



{i^\{r)Mr')) - ^(V^I(r)V^i(r')) + x(V^|(r)Vx(r')) 



^(^Ur')V|(r)) + ^(^|(r)V^i(r')). 

(B5) 



In this way, we can rewrite the last lines of Eq. (B3 1 as: 



lim(V.-Vw)[i(G^'(r,r')) 

am r^r' \ / 21 



-1 G^(r,r') -1 G^(r,r') 



34 



.(G^(r,r')) 



43 



= lim (V^ - Vr')Tr{p3 

am r^r' 

xdiag(T2,T2*) X G^(r,r')}. 



(B6) 



For the x and z-components, one replaces T2 with ti and rs, 
respectively. 



We here show how the matrix structure in the analytical ex- 
pression Eq. (15 I for the spin-current is obtained in the qua- 
siclassical approximation. The starting point is the quantum 



APPENDIX C: PAULI MATRICES 



The Pauli matrices that are used in this paper are 
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